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In  this  paper,  a  thermoelectric  generator  (TEG)  energy  harvesting  system  with  a  temperature- 
sensor-based  maximum  power  point  tracking  (MPPT)  method  is  presented.  Conventional  MPPT 
algorithms  for  photovoltaic  cells  may  not  be  suitable  for  thermoelectric  power  generation  because  a 
significant  amount  of  time  is  required  for  TEG  systems  to  reach  a  steady  state.  Moreover,  complexity 
and  additional  power  consumption  in  conventional  circuits  and  periodic  disconnection  of  power  source 
are  not  desirable  for  low-power  energy  harvesting  applications.  The  proposed  system  can  track  the  vary¬ 
ing  maximum  power  point  (MPP)  with  a  simple  and  inexpensive  temperature-sensor-based  circuit  with¬ 
out  instantaneous  power  measurement  or  TEG  disconnection.  This  system  uses  TEG’s  open  circuit  voltage 
(OCV)  characteristic  with  respect  to  temperature  gradient  to  generate  a  proper  reference  voltage  signal, 
i.e.,  half  of  the  TEG’s  OCV.  The  power  converter  controller  maintains  the  TEG  output  voltage  at  the  refer¬ 
ence  level  so  that  the  maximum  power  can  be  extracted  for  the  given  temperature  condition.  This  feed¬ 
forward  MPPT  scheme  is  inherently  stable  and  can  be  implemented  without  any  complex  microcontroller 
circuit.  The  proposed  system  has  been  validated  analytically  and  experimentally,  and  shows  a  maximum 
power  tracking  error  of  1.15%. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Recently,  investigation  into  sustainable  alternative  energy 
sources  has  been  significantly  expanded.  In  addition  to  the  large- 
scale  renewable  energy  systems  using  wind  turbines  and  photovol¬ 
taic  (PV)  panels,  the  harvesting  of  small  energy  from  various  free 
energy  sources,  such  as  vibrations  from  a  bridge,  waste  heat  from 
internal  combustion  engines  and  brick  kilns,  and  electricity  from 
biological  molecules,  has  been  investigated  in  regard  to  direct 
power  generation  for  low-power  loads  or  system  efficiency 
improvement  through  energy  recovery  [1,2].  Among  these 
approaches,  heat  energy  harvesting  via  thermoelectric  generators 
(TEGs)  converts  thermal  energy  to  electrical  energy  without  any 
mechanical  moving  part.  It  is  small,  light,  reliable,  eco-friendly, 
and  promising  for  applications  such  as  automotive  waste  heat 
recovery  and  wireless  sensor  network  power  supply  due  to  its  long 
lifetime  and  high  reliability  [3-6].  Although  it  has  suffered  from 
low  efficiency  and  high  cost,  the  feasibility  of  TEG  as  an  energy 
source  is  being  improved  by  recent  research  progress  in  areas  such 
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as  material  efficiency  [7-9]  and  energy  harvesting  system  develop¬ 
ment  [10,11]. 

While  most  research  efforts  have  focused  on  enhancing  thermo¬ 
electric  material  properties,  power  conditioning  circuits  optimized 
for  TEG  systems  have  not  been  investigated  extensively.  Many 
studies  have  suggested  optimal  electrical  load  matching  conditions 
for  maximized  power  generation,  and  utilized  maximum  power 
point  tracking  (MPPT)  techniques  originally  developed  for  PV 
systems,  such  as  Perturb  and  Observe  (P&O)  and  Incremental 
Conductance  (IC)  [12-14].  Although  maximizing  the  energy  har¬ 
vest  from  TEGs  by  means  of  an  MPPT  technique  is  essential  for 
TEG  applications  to  compensate  the  TEG’s  low  conversion  effi¬ 
ciency,  conventional  MPPT  techniques  for  PV  systems  may  not  be 
optimal  for  TEG  systems.  Unlike  a  PV  system,  a  TEG  system 
requires  a  certain  amount  of  time  to  reach  a  steady-state;  hence, 
fast  perturbations  on  output  power  in  P&O  and  IC  algorithms  could 
lead  to  off-MPP  operations  due  to  the  inaccurate  instantaneous 
TEG  power  measurement.  Moreover,  these  hill  climbing  techniques 
rely  on  fast  measurements  and  complex  computation,  which  in 
turn  require  additional  hardware  and  software  that  consume  mea¬ 
surable  amount  of  power.  The  fractional  voltage  technique  that 
samples  the  open  circuit  voltage  (OCV)  of  the  TEG  while  in  opera¬ 
tion  is  also  widely  used  [15-17].  But,  it  also  suffers  from  the  large 
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system  time  constant  that  is  required  to  sample  precise  OCV  for 
the  given  temperature  differential.  Furthermore,  a  microcontroller 
and  analog  switch  are  required  to  disconnect  the  TEG  from  the 
circuit  and  sample  OCV,  which  can  be  costly  and  troublesome 
especially  in  low-power  TEG  systems. 

In  this  paper,  a  temperature-sensor-based  MPPT  technique  is 
proposed  for  TEG  energy  harvesting  systems.  The  proposed  system 
can  track  variations  of  MPP  that  are  due  to  temperature  changes 
with  a  simple  and  inexpensive  temperature-sensor-based  circuitry 
that  can  avoid  the  disadvantages  of  conventional  techniques,  such 
as  current  measurement,  source  disconnection,  and  power  pertur¬ 
bation.  Instead,  TEG’s  OCV  characteristic  with  respect  to  tempera¬ 
ture  difference  between  source  and  environment  is  built  into  the 
circuit;  then  a  proper  command  voltage  signal  -  i.e.,  half  of  the 
TEG’s  OCV  for  the  given  temperature  differential  -  is  generated. 
The  command  is  used  by  a  hysteresis  controller  to  generate  a  gate 
signal  for  the  power  converter  to  harvest  power  from  the  TEG  at 
the  MPP  of  the  given  temperature  condition  without  perturbation 
or  interruption.  The  proposed  MPPT  circuitry  can  be  implemented 
with  a  minimal  number  of  discrete  parts  and  without  any  complex 
microcontroller  system.  Moreover,  the  system  is  inherently  stable 
because  it  operates  based  on  the  feedforward  model  without  any 
instantaneous  power  feedback  loop.  The  proposed  scheme  has 
been  validated  analytically  and  experimentally,  and  demonstrated 
successful  performance. 

2.  TEG  energy  harvesting 

2.2.  Thermal  model  of  thermoelectric  module 


temperatures  of  the  heat  source,  hot  surface,  cold  surface,  and 
ambient  air,  respectively.  The  energy  conservation  equations  can 
be  further  expanded  as 

0.H  =  K(Th  -  Tc)  +  OcItecTh  - 1 12TECRint  (3) 

Qc  =K(Th-  Tc)  +  ocIjecTc  +  \  TTECRi nt  (4) 

where  I<  is  the  thermal  conductance  of  TEG  material,  a  is  the 
Seebeck  coefficient  of  TEG  material,  ITEc  is  the  TEG  output  current, 
and  Rint  is  the  internal  electrical  resistance  of  the  TEG  material 
[18].  The  Seebeck  coefficient  and  the  internal  resistance  have  been 
shown  to  be  nearly  constant  over  the  operating  range  of  the  TEG 
[15].  The  output  power  of  the  TEG  module  will  be  given  as 

P  =  Qh~Qc  (5) 

Because  Qpeitier  and  Qjh  are  functions  of  the  TEG  output  current 
as  can  be  seen  in  (3)  and  (4),  the  TEG  output  power  will  vary  as  the 
output  current  changes  even  with  a  constant  temperature  differ¬ 
ence  between  TH  and  Tc.  Furthermore,  TH  and  Tc  will  also  change 
with  the  amount  of  output  current  for  the  given  source  and  ambi¬ 
ent  temperatures  Ts  and  T^,  because  QH  and  Qc  are  functions  of  the 
output  current.  So  far,  many  thermoelectric  system  analyses  have 
assumed  that  change  of  temperatures  by  electrical  current  is  neg¬ 
ligible,  and  have  used  a  conventional  maximum  power  output  con¬ 
dition  ( Rint  =  Rext).  Several  researchers  addressed  the  inaccuracy  of 
this  approach  in  real  applications  and  suggested  new  optimum 
operation  conditions  18-21]. 


Thermoelectric  modules  directly  generate  electricity  propor¬ 
tional  to  the  temperature  difference  applied  across  the  legs.  A  typ¬ 
ical  schematic  of  a  TEG  system  and  its  thermal  circuit  model  is 
shown  in  Fig.  1.  The  heat  flow  in  the  system  can  be  represented 
by  the  following  two  energy  conservation  equations. 


Qh 


Qcond  +  Qpeitier  QjH 


Ts-Th 

Vh 


a) 


Qc  —  Qcond  +  Qpeitier  +  QjH  — 


Tc  ~  T^ 

^C 


(2) 


where  QH  and  Qc  denote  the  amount  of  heat  transfer  at  the  hot  and 
cold  surfaces  of  the  thermoelectric  module,  QCond ,  Qpeitier ,  and  QJH  are 
heat  conduction,  Peltier  heating/cooling,  and  Joule  heating  inside 
the  thermoelectric  material,  respectively.  \j/H  and  i/rc  are  the  thermal 
resistance  of  the  hot  and  cold  sides.  Ts,  TH,  Tc,  and  are  the 
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Fig.  1.  Schematic  and  thermal  circuit  model  of  a  typical  thermoelectric  module. 


2.2.  Electrical  model  and  maximum  power  point 

The  electrical  power  converted  from  the  thermal  power  by  TEG 
will  be  given  as 

P  =  V  TEG  I  TEG  (6) 

where  VTEC  and  ITEG  are  the  voltage  and  current  at  TEG’s  output  ter¬ 
minals,  respectively.  Taking  into  account  the  thermal  factors 
described  in  Section  2.1,  the  TEG  can  be  modeled  as  a  dependent 
voltage  source  Vs  with  an  internal  resistance  Rint,  as  shown  in 
Fig.  2.  The  external  resistance  Rext  represents  the  electrical  load  that 
draws  power  from  TEG. 

The  TEG  will  operate  at  a  point  on  the  V-I  curve  or  power  curve. 
It  can  also  be  associated  with  the  load  resistance  seen  by  the  TEG. 
In  the  conventional  system  model  where  the  internal  voltage  Vs  is 
constant  and  independent  of  output  current,  it  can  be  analytically 
derived  that  the  MPP  occurs  when  Rint  =  Rext  [15].  However,  it  has 
been  shown  that  this  impedance-matched  MPP  condition  is  valid 
only  when  the  temperatures  at  both  ends  of  the  thermoelectric 
legs  are  kept  constant  [22,23].  This  assumption  cannot  hold  in 
practical  applications,  because  temperatures  at  both  ends  of  the 
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thermoelectric  legs  that  determine  the  internal  source  voltage  in 
the  equivalent  circuit  are  affected  by  multiple  factors,  including 
conduction  heat  transfer  through  heat  sink,  energy  transfer  by 
Joule  heating,  the  Peltier  effect,  and  thermal  conduction  in  the 
TEG  module.  Even  for  a  constant  temperature  differential,  the 
internal  voltage  will  change  as  load  current  varies,  which  makes 
the  MPP  deviate  from  where  the  impedances  are  matched.  For 
example,  if  Vs  decreases  as  output  current  increases  (Vs2),  the  V-I 
curve  and  power  curve  change  from  the  constant  V*  case  (Vsi) 
and  its  MPP  (MPP2)  will  shift  to  where  Rextf  >  Rint,  as  shown  in 
Fig.  3.  However,  it  can  be  seen  from  the  V-I  lines  that  the  MPPs 
are  still  at  OCV/2. 

New  forms  of  load  matching  conditions  have  been  suggested 
[18,19],  but  such  conditions  can  only  be  applied  when  system 
parameters,  such  as  heat  sink  or  geometry  of  TEG,  are  carefully 
chosen.  Moreover,  load  conductance  is  difficult  to  match  when  a 
non-constant  electric  load  is  directly  connected;  however,  because 
the  V-I  curve  is  near  linear  for  most  applications  [24,25],  especially 
when  source  temperature  is  kept  constant  [26],  the  MPP  will  be  at 
half  of  the  OCV  regardless  of  RextlRint  relationships. 


techniques,  this  perturbation-based  approach  is  not  suitable  for 
TEG  applications  due  to  TEG’s  slow  response  to  the  power 
perturbation. 

The  virtual  conductance  measurement  technique  uses  voltage 
and  current  measurements  to  calculate  instantaneous  load  con¬ 
ductance  {G0  =  ItegIVjeg)  and  compare  it  to  the  pre-determined 
internal  conductance  (G/=  1/Pint)  to  find  the  MPP  [34].  The  power 
converter  controls  the  TEG  output  voltage  VtEG  to  operate  at  the 
load  matching  condition;  however,  as  shown  in  Fig.  3,  this 
approach  cannot  be  applied  to  the  system  that  the  internal  voltage 
source  changes  as  output  current  varies  because  the  load  matching 
is  not  the  MPP  condition. 

The  fractional  OCV  method  uses  the  relationship  between  MPP 
voltage  and  OCV  [16,17,35].  The  OCV  is  sampled  periodically  to 
determine  the  MPP  voltage  for  the  given  condition;  however,  the 
power  source  should  be  electrically  disconnected  by  a  switch  or 
by  shutting  off  the  power  converter  for  the  OCV  measurement 
while  in  operation.  Hence,  temporary  power  loss  and  complex 
implementation  are  drawbacks  if  this  approach.  For  TEG  applica¬ 
tions,  the  large  time  constant  of  the  temperature  change  would 
also  be  problematic  for  precise  OCV  measurement. 


2.3.  Conventional  MPPT  techniques 

The  Perturb  and  Observe  (P&O)  technique  is  widely  used  in  TEG 
applications  as  well  [27-29].  It  samples  the  instantaneous  power 
production  of  the  generator  after  a  perturbation  at  discrete  inter¬ 
vals  on  either  side  of  the  current  operating  point,  and  moves  the 
operating  point  toward  that  of  higher  power;  however,  the  P&O 
technique  responds  poorly  to  a  fast-changing  power  output  [30] 
and  it  is  inefficient  at  the  steady  state  due  to  the  fact  that  the  con¬ 
stant  perturbation/sampling  and  high-frequency  fluctuation  of  the 
operating  point  cause  a  harmonic  around  the  MPP  [31]. 

Incremental  conductance  (IC)  is  a  method  that  can  avoid  the 
oscillations  at  MPP  in  P&O  systems  [32].  The  IC  method  compares 
the  incremental  conductance  {Gd  =  dl/dV)  of  the  generator  against 
the  negative  value  of  instantaneous  conductance  (Gs  =  -J/V).  At 
any  point  on  the  power  curve  before  the  MPP,  Gd  is  greater  than 
Gs  and  the  operational  point  is  moved  toward  the  higher  voltage. 
At  an  operating  point  after  the  MPP,  Gs  is  greater  than  Gd  and  the 
operational  point  is  moved  toward  a  lower  voltage.  The  MPP  is 
reached  when  Gd  =  Gs,  and  IC  does  not  search  unless  there  is  a 
change  in  current.  The  IC  algorithm  shows  better  performance  in 
terms  of  noise  and  system  dynamics  compared  to  P&O,  but  worse 
in  rapidly  changing  source  conditions,  and  there  would  be  a 
stability  issue  with  a  high  perturbation  rate  [33].  Like  the  P&O 


3.  Proposed  MPPT  system 

In  this  paper,  a  model-based  MPPT  controller  that  estimates 
OCV  using  two  temperature  measurements  is  developed  to  avoid 
the  aforementioned  disadvantages  of  conventional  techniques. 
The  proposed  system  determines  the  operating  point  at  MPP, 
which  is  half  of  the  OCV,  in  a  feedforward  manner  using  tempera¬ 
ture  measurements  without  any  feedback  information  of  actual 
power  or  output  resistance.  This  results  in  continuous  power 
harvesting  at  MPP  that  does  not  require  any  perturbation  or  inter¬ 
ruption.  Furthermore,  complex  hardware  and  software  require¬ 
ments  and  complicated  thermal-electrical  considerations  are  not 
necessary,  because  the  MPP  is  at  half  of  the  OCV  for  most  of  the 
applications  due  to  the  linear  V-I  characteristics. 


3.1.  OCV  estimation  using  temperature  measurements 


From  the  thermal  model  shown  in  Fig.  1,  temperature  drop 
across  thermoelectric  legs  in  the  open  circuit  condition  when 
Imfc  =  0  can  be  obtained  as  follows. 


(Th  -  Tc)oc  = 


yi< 


'Fh  +  'Fc  +  'L/K 


(Ts  -  TJ  = 


(Ts-y 


KWh+KVc  +  I 


(7) 


Fig.  3.  Comparison  of  operations  with  constant  (Vsi)  and  output-current-dependent  (Vs2)  internal  voltage  source  for  a  given  temperature  difference  between  both  ends  of 
thermoelectric  legs.  It  can  be  seen  that  the  MPP2  shifts  to  where  Rext  ^  Rint,  but  still  occurs  at  OCV/2.  (a)  TEG  output  currents  and  power  vs.  TEG  output  voltage,  (b)  Output 
power  vs.  Rext/Rint  ratio. 
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The  OCV  is  directly  proportional  to  the  temperature  difference 
at  the  open  circuit  condition,  and  the  proportional  coefficient  will 
be  the  Seebeck  coefficient  5  of  a  TEG  module. 

Voc-sir.-rcW-i^iL  (s, 

The  thermal  resistances  depend  on  contact  resistances,  pad 
resistance,  and  heat  sinks,  among  others,  but  these  values  are  quite 
difficult  for  end  users  to  evaluate  and  must  be  identified  through 
calibration.  The  Seebeck  coefficient  S  of  a  TEG  module  that  contains 
N  p-  and  n-type  legs  will  be  given  as  Not ,  where  a  is  the  material 
Seebeck  coefficient. 

It  can  be  seen  from  (7)  that  the  OCV  can  be  represented  as  a  lin¬ 
ear  function  of  external  temperatures  of  Ts  and  T^. 

Voc  =  K(TS  -  Too)  (9) 

where  k  is  the  proportional  constant,  which  depends  on  thermal 
resistances  and  module  properties.  Although  it  is  difficult  to  analyt¬ 
ically  calculate  k  because  module  properties  are  not  usually  pro¬ 
vided  by  TEG  manufacturers,  the  constant  k  can  be  readily 
determined  experimentally  by  measuring  OCV  for  the  given  tem¬ 
perature  differences  between  Ts  and  T^.  It  does  not  depend  on 
the  output  current  because  it  is  measured  in  no-load  conditions. 
Furthermore,  Ts  and  I ^  are  easily  accessible  unlike  TH  and  Tc. 

3.2.  Feedforward  MPPT  controller  for  uninterrupted  energy  harvesting 

The  proposed  MPPT  controller,  i.e.,  a  TEG  model/MPP  reference 
generator,  implements  an  energy  harvesting  system  in  conjunction 
with  a  hysteresis  controller  to  generate  a  power  converter  gate  sig¬ 
nal  and  a  boost  converter  for  energy  capture.  The  proposed  energy 
harvesting  system  extracts  and  stores  the  maximum  attainable 
energy  for  the  given  operating  conditions.  A  schematic  of  the  over¬ 
all  system  is  shown  in  Fig.  4. 

TEG’s  OCV  is  linear  in  terms  of  temperature  differential  as  can 
be  seen  in  (9),  and  it  has  been  shown  that  MPP  is  at  half  of  the 
OCV.  The  reference  voltage  for  MPP  will  therefore  be 

VREF=^(TS-T00)=^AT  (10) 


This  TEG  model  -  i.e.,  OCV  vs.  temperature  difference  -  can  easily  be 
modeled  with  a  differential  amplifier  and  temperature  sensors  as 
shown  in  Fig.  4,  which  enables  a  feedforward  tracking  of  the  MPP 
based  on  instantaneous  temperature  measurements.  The  output 
voltage  of  the  model  circuit  is  given  as  follows: 

Vrbf=^AVt  (11) 


where  AVT  is  the  voltage  difference  of  temperature  sensor  output 
for  Ts  and  T^  (i.e.,  VHOt~  VCOld)-  Using  a  circuit  with  potentiome¬ 
ters,  the  OCV  characteristic  of  any  TEG  module  can  be  modeled, 
and  the  OCV/2  for  any  temperature  differential  can  be  readily  gen¬ 
erated  for  the  power  conversion  stage.  The  temperature  measure¬ 
ment  circuit  is  much  easier  to  implement  and  less  expensive  than 
the  ones  used  for  current  measurement  or  TEG  disconnection  which 
require  a  microcontroller  system  for  A/D  conversion  and  MPPT  algo¬ 
rithm  implementation,  which  is  a  significant  advantage  especially 
for  low-power  systems.  Continuous  power  harvesting  without  any 
perturbation  or  interruption  would  also  be  a  substantial  advantage 
of  the  proposed  scheme  for  TEG  applications  that  have  a  large  time 
constant  for  temperature  changes.  Moreover,  it  can  be  implemented 
with  an  OP  amp  and  four  resistors,  which  is  much  simpler  than  the 
conventional  microcontroller-based  MPP  searching  algorithms. 

To  maintain  the  TEG  output  voltage  at  the  commanded  MPP,  the 
hysteresis-controller-based  gate  drive  circuitry  is  used  [36,37].  The 
hysteresis  control  circuit  generates  an  appropriate  gating  signal  so 
that  the  power  converter  can  maintain  the  TEG  output  voltage  VTEG 
at  the  reference  voltage  of  the  given  temperature  condition,  which 
tracks  the  varying  MPP.  A  hysteresis  controller  can  easily  be  imple¬ 
mented  with  a  comparator  and  two  resistors,  as  shown  in  Fig.  4. 

The  hysteresis  controller  has  two  thresholds  {VThH  and  VThL) 
generated  by  VREF  to  compare  with  the  input  voltage  VTEG. 


VThH  = 


R3+R4 

R4 


Vref~mv°l 


(12) 


VjhL  = 


R3  +  R4 
R4 


Vref  — 


R3 

R4VoH 


(13) 


where  VREF  is  the  output  voltage  of  the  TEG  model  circuit,  which  is 
OCV/2  for  the  given  temperature  difference,  and  VoH  and  VoL  are  the 


Energy  harvesting  boost  converter 


Fig.  4.  Overall  schematic  of  the  proposed  energy  harvesting  system.  TEG  model  feeds  the  hysteresis  controller  with  an  MPP  voltage  reference.  Hysteresis  controller  controls 
boost  converter  and  maintains  energy  harvesting  at  MPP. 
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high  and  low  voltage  of  the  comparator’s  binary  output  VG,  respec¬ 
tively.  When  the  TEG  output  voltage  reaches  the  upper  threshold 
VThH,  the  hysteresis  controller  turns  on  the  main  power  MOSFET  Q. 
in  the  boost  converter,  (VG  =  VoH ,  logic  high),  to  start  extracting 
power  from  the  TEG.  In  this  mode,  VTEC  decreases  as  output  current 
increases  and  the  energy  is  stored  in  the  inductor  of  the  boost  con¬ 
verter.  When  the  VTEC  hits  the  lower  threshold  VThL,  the  hysteresis 
controller  turns  off  Q.  (VG  =  VoL,  logic  low)  and  the  extraction  stops. 
Then,  energy  stored  in  the  inductor  is  discharged  to  the  output 
capacitor.  In  this  mode,  VTEC  increases  and  when  it  is  recovered  to 
VThH,  the  extraction  resumes.  The  output  of  the  hysteresis  controller 
can  be  seen  in  Fig.  5.  The  high  and  low  logic  levels  of  VG  to  turn  on 
and  off  the  power  switch  can  vary  by  circuit  configuration. 

If  we  pick  an  R4  value  significantly  larger  than  R3  ( R4  »  R3),  the 
thresholds  will  be  given  as  follows  for  a  general  comparator  with  a 
dual  power  supply  {—Vsat  <  Vo  Vsatl 


DO 

VThH  =  VREF  +  M  Vsat 

(14) 

DO 

VThL  =  VREF-—Vsat 

(15) 

The  control  circuit  can  be  much  simpler  if  a  single  supply  compar¬ 
ator  (0  <  V0  <  Vsat )  is  used.  Then,  the  upper  threshold  becomes 

Vthh  -  VREE  (16) 

The  gate  signal  VG  changes  between  zero  and  VSAT  if  the  circuit  is 
implemented  with  a  single  supply  op  amp. 

By  repeating  these  two  operational  modes,  VTEG  can  be  confined 
in  the  band  between  VThH  and  VThL  around  MPP.  The  periods  of 
these  two  modes,  as  well  as  the  switching  frequency,  are  deter¬ 
mined  by  the  TEG  output  power  and  power  converter  parameters. 
The  operation  of  the  TEG  with  a  hysteresis  controller  is  shown  in 
Fig.  6. 

4.  Results  and  discussion 

4A.  Experimental  setup 

A  diode-based  boost  converter  is  selected  for  the  energy  har¬ 
vesting  experiment  in  this  paper  due  to  its  simplicity,  but  any  type 
of  power  converter  can  be  integrated  with  the  proposed  control 
scheme.  A  low-power  single-supply  operational  amplifier 
(LM392N,  Texas  Instruments,  TX),  a  low-power  comparator 
(MCP6542,  Microchip  Technology,  AZ),  and  generic  resistors  and 
potentiometers  have  been  used  to  implement  feedforward  TEG 
model  and  a  hysteresis  controller,  respectively.  For  the  boost  con¬ 
verter,  a  450pH  inductor  (RC-7,  Triad  Magnetics,  CA),  an  N-channel 
MOSFET  (NTD4906N-35G,  ON  Semiconductor,  AZ),  a  Schottky 
diode  (BAT46,  STMicroelectronics,  Switzerland),  and  a  3F,  2.5  V 
supercapacitor  (M1020-2R5305-R,  Cooper  Bussmann,  MO)  were 
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Fig.  5.  Thresholds  VThH,  VThH  and  output  VG  of  hysteresis  controller. 


used.  A  generic  potentiometer  (1  MQ)  was  used  to  control  the  load 
current.  The  implemented  MPPT  control  circuit  and  power  con¬ 
verter  can  be  seen  in  Fig.  7. 

The  hot  side  of  the  TEG  module  was  mounted  on  an  aluminum 
block  on  a  temperature-controlled  hotplate  (HP131225Q,  Thermo 
Scientific,  MA),  and  the  cold  side  was  attached  to  a  generic  CPU 
cooling  heat  sink/fan  assembly.  A  fire  brick  was  inserted  between 
the  hotplate  and  the  heat  sink  to  minimize  the  effect  of  convection 
heat.  For  temperature  measurement,  temperature  sensors 
(MCP9700,  Microchip  Technology,  AZ)  were  used.  MCP9700  has  a 
linear  temperature  slope  (10  m V/°C)  over  the  whole  operating 
range  unlike  thermistors  and  significantly  higher  output  voltage 
(1.75  V  at  125  °C)  than  thermocouples.  Furthermore,  this  sensor 
has  good  measurement  accuracy  (±4  °C  max)  and  draws  low  cur¬ 
rent  (6  pA  typical),  so  it  can  be  readily  used  for  low-tempera- 
ture-difference  and  low-power  TEG  applications.  Thermal 
compound  (RG-ICFN-200G-B1,  Cooler  Master,  CA)  and  thermal 
adhesive  (ASTA-7G,  Arctic  Silver,  CA)  were  used  for  the  TEG  and 
temperature  sensor  mount,  respectively.  A  current  meter  (AEMC 
K110,  Chauvin  Arnoux,  NH)  with  a  lOmV/mA  gain  was  used  to 
measure  the  TEG  output  current.  Two  thermocouples  and  multim¬ 
eters  were  used  to  measure  temperatures;  these  measurements 
were  for  reference  and  were  not  used  for  control.  For  the 
experiments  in  this  paper,  a  Bi-Te  based  TEG  module  (G2-30- 
0313,  Tellurex,  MI)  was  used.  The  experimental  setup  is  shown 
in  Fig.  8. 

4.2.  Results 

To  determine  the  TEG  model  circuit  parameters,  the  experimen¬ 
tal  TEG  system’s  OCV  in  terms  of  the  temperature  difference  AT 


Fig.  7.  Energy  harvesting  system:  MPPT  controller,  hysteresis  controller,  and  power 
converter.  Load  resistor  and  ambient  temperature  sensor  are  also  shown. 
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Fig.  8.  Experimental  setup. 


between  Ts  and  T ^  was  characterized,  and  the  result  is  shown  in 
Fig.  9.  A  linear  characteristic  can  be  seen  as  expected  in  (9),  and 
the  coefficient  k  was  found  to  be  0.0154  [V/°C].  The  OCV  character¬ 
istic  with  respect  to  the  temperature  sensor  output  was  also 
measured,  and  it  is  characterized  as 

Voc  =  0.7  4AVt  (17) 

where  AVT  is  the  output  voltage  difference  of  temperature  sensors 
for  Ts  and  T^.  Hence,  the  command  voltage  for  the  MPP  of  the  given 
voltage  differential  (i.e.,  temperature  differential)  will  be  given  as 
follows: 


ambient  temperature  change  was  negligible  during  testing.  Around 
15  data  points  were  taken  in  load  current  range  from  0  to  600  mA, 
and  enough  time  was  given  to  stabilize  the  output.  Data  points  and 
linear  regression  lines  can  be  seen  in  Fig.  10.  The  linear  expressions 
for  each  V-I  characteristic  are  shown  in  Table  1,  and  OCVs  esti¬ 
mated  from  V-I  characteristics  are  compared  to  experimentally 
measured  ones  to  verify  the  expressions. 

The  proposed  MPPT  scheme  was  tested  with  the  same  temper¬ 
ature  conditions  that  were  used  for  V-I  characterization.  When  the 
temperature  condition  is  changed,  the  TEG  model  is  supposed  to 
feed  the  appropriate  command  voltage  (i.e.,  MPP  voltage  for  the 
new  temperature  differential)  to  the  hysteresis  controller  so  that 
the  power  converter  can  always  extract  the  maximum  power  in 
any  temperature  condition.  The  performance  of  the  proposed  sys¬ 
tem  has  been  verified  by  monitoring  instantaneous  waveforms  of 
TEG  output  voltage  and  current,  the  output  signal  of  the  hysteresis 
controller,  and  the  power  converter  output  voltage  using  a  digital 
oscilloscope.  A  screen  capture  for  the  operation  when  AT  =  17  °C 
is  shown  in  Fig.  11.  It  can  be  seen  that  the  TEG  is  operating  at 
157  mV  (CHI)  and  70.5  mA  (CH3,  ImA/lOmV),  which  is  the 
MPP  for  the  given  temperature  condition.  The  operations  at  other 
temperature  conditions  were  measured  in  the  same  manner,  and 
the  result  of  MPPT  operation  is  shown  along  with  the  voltage- 
power  curves  generated  from  V-I  characteristics  in  Fig.  12.  The 
solid  dots  denote  the  operating  points  controlled  by  the  proposed 
controller  at  each  given  temperature  condition,  and  comparisons 
to  the  calculated  operating  conditions  from  the  V-I  expressions 
are  shown  in  Table  1.  Clearly,  the  proposed  controller  tracks  the 
varying  MPP  fairly  well. 

The  hysteresis  voltage  band  was  set  to  100  mV  using  a  fixed 
30  Q  resistor  and  3.25  MQ  potentiometers  for  R3  and  R4,  respec¬ 
tively.  It  should  be  noted  that  the  band  could  be  larger  than  the 
theoretical  value  because  of  the  propagation  delay  of  the 
comparator. 


4.3.  Discussion 


Vcmd  =  0.5  xVoc  =  0.37 AVt  (18) 

The  command  voltage  is  also  shown  in  Fig.  9.  This  TEG’s  MPP 
characteristic  in  terms  of  the  sensor  output  voltage  difference  is 
implemented  in  the  model  circuit  using  20  kQ  and  7.42  1<Q  for  R i 
and  R2,  respectively.  These  two  resistances  are  the  only  parameters 
to  be  determined  for  TEG  model.  Although  any  resistance  values  can 
be  selected  as  long  as  the  ratio  is  maintained,  higher  values  would 
be  preferable  due  to  low  power  consumption. 

The  V-/  characteristic  of  the  overall  system  is  identified  for  dif¬ 
ferent  load  and  temperature  conditions.  Five  different  temperature 
conditions  were  given:  AT  =  17  °C,  45  °C,  79  °C,  95  °C  and  122  °C. 
The  source  temperature  was  regulated  by  the  hotplate,  and 


The  MPPT  circuit  suggested  in  this  paper  estimates  OCV  from 
source  and  ambient  temperatures  that  can  be  easily  measured 
and  extracts  energy  at  half  of  the  estimated  OCV.  One  of  the  advan¬ 
tages  of  the  proposed  technique  over  other  MPPT  algorithms  is  that 
it  can  be  implemented  with  much  less  energy  consumption  and 
complexity.  As  the  power  and  voltage  generated  from  a  TEG  mod¬ 
ule  are  quite  small,  even  a  small  voltage  drop  and  power  loss  across 
the  shunt  resistor  can  be  detrimental  and  may  not  be  economically 
justified.  The  same  is  true  of  any  additional  power  consumption  by 
complex  power  management  circuit  using  A/D  converter  and 
microcontroller.  Although  the  MPPT  systems  based  on  the  conven¬ 
tional  fractional  voltage  is  fast  and  reliable,  the  technical  challenge 
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Fig.  9.  Experimental  OCV  measurement  with  respect  to  temperature  difference  and 
temperature  sensor  output  voltage  difference  between  Ts  and  T^.  Command  voltage 
for  power  converter  controller,  OCV/2,  is  also  shown. 


Fig.  10.  V-I  characteristics  of  the  TEG  used  in  this  paper  with  different  load  and 
temperature  conditions. 
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Table  1 

Experimental  result:  (1)  V-I  characterizations.  (2)  OCV  measurements.  V0c '■  estimated  OCV  from  V-I  curve  expression,  V0c ’  experimentally  measured  OCV.  (3)  Operating  point 
comparisons.  VTEG,  I  teg  and  P:  calculated  output  current  and  power  at  MPP  from  V-I  characteristics.  VTEG ,  I  teg  and  P:  actual  TEG  output  voltage,  current,  and  power  controlled  by 
the  proposed  system.  (4)  Performance.  VERR  and  PERR  denote  the  error  between  VTEG  and  VTEG ,  PTEG  and  PTEG. 


A T(°C) 

System  characterization 

V-I  curve  expression  (mA,  mV) 

OCV  measurements 

Calculated  operating  point 

Actual  operating  point 

Performance 

V0C  (mV) 

V0c  (mV) 

Err  (%) 

Vteg  (mV) 

I teg  (mA) 

P  (mW) 

Vteg  (mV) 

Iteg  (mA) 

P  (mW) 

VERR  (%) 

Perr  (%) 

17 

i  =  -0.4547  z/+ 141.83 

312.0 

311.0 

0.32 

155.0 

70.9 

11.1 

157.0 

70.5 

11.1 

-1.29 

0.08 

45 

i  =  -0.4587^+339.78 

740.8 

735.0 

0.78 

367.5 

169.9 

62.9 

370.0 

172.0 

63.6 

-0.68 

1.15 

79 

i  =  -0.4301  v+  523.05 

1216.2 

1210.0 

0.51 

605.0 

261.5 

159.0 

615.0 

258.5 

159.0 

-1.65 

-0.03 

95 

i  =  -0.4505  v+  629.95 

1534.6 

1530.0 

0.30 

765.0 

315.0 

241.7 

769.0 

316.0 

243.0 

-0.52 

0.55 

122 

i  =  -0.3661  722.39 

1973.2 

1960.0 

0.67 

986.6 

361.2 

356.4 

986.6 

365.0 

358.5 

0.36 

0.68 
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Fig.  11.  Instantaneous  TEG  voltage,  current,  power  converter  gate  signal,  power 
converter.  TEG  current  (CH3)  is  measured  with  1  mA/10  mV  current  probe.  The  TEG 
is  operating  at  157  mV  (CHI)  and  70.5  mA(CH3, 1  mA/10  mV),  which  is  the  MPP  for 
the  given  temperature  condition. 


Fig.  12.  Experimental  result  of  MPPT  control.  Solid  dots  denote  the  controlled 
operating  points  by  the  proposed  system  at  each  given  temperature  condition. 

lies  in  how  to  determine  OCV  efficiently.  Commercially  available 
OCV-based  MPPT  circuits  periodically  cut  electrical  current  from 
TEGs  and  sample  OCV,  which  has  critical  issues  with  TEGs.  If  cur¬ 
rent  is  interrupted  for  a  short  period  in  order  to  achieve  a  faster 
tracking  performance,  it  may  not  be  long  enough  for  TEGs  to  reach 
thermal  equilibrium  and  correct  OCV.  Similarly,  P&O  and  IC  meth¬ 
ods  have  the  same  issue  as  the  temperature  cannot  respond 
quickly  enough  to  the  perturbations.  The  change  in  power  after 
perturbation  can  be  instantaneously  larger  than  its  steady-state 
value.  It  could  even  take  a  few  minutes  for  some  systems  to  reach 
a  steady-state,  and  such  long  disconnection  is  not  desirable  for 
energy  harvesting  applications. 

The  proposed  feedforward  model  and  the  hysteresis  controller 
can  be  implemented  with  a  minimal  number  of  discrete  parts, 


and  the  system  can  be  further  simplified  if  the  amplifier  and  com¬ 
parator  are  built  into  a  single  chip.  The  proposed  scheme  does  not 
require  any  current  measurement  or  complex  microcontroller 
system,  which  makes  the  proposed  scheme  very  straightforward, 
cost-effective,  and  reliable.  Although  only  an  energy  harvesting 
converter  that  captures  the  maximum  available  energy  from  TEG 
was  built  in  this  study,  a  secondary  power  converter  can  easily 
be  implemented  for  the  additional  output  voltage  control,  such 
as  constant  voltage  output  or  battery  charging,  to  provide  more 
flexible  operation. 

The  suggested  circuit  can  steadily  find  an  MPP  even  when  the 
energy  harvester  is  exposed  to  unstable  source/ambient  tempera¬ 
ture,  such  as  waste  heat  recovery  from  a  cook  stove  or  a  wearable 
energy  harvester,  as  well  as  under  controlled  cooling  conditions. 
One  should  note  that  the  proposed  circuit  responds  to  the  change 
in  source/ambient  temperatures  and  cannot  track  the  MPP  varia¬ 
tion  due  to  the  change  of  k,  the  heat  dissipation  characteristic  of 
the  energy  harvesting  system;  however,  although  the  convective 
heat  transfer  coefficient  is  a  function  of  temperature,  the  overall 
change  in  k  is  not  significant  for  stationary  applications.  A  new 
MPPT  strategy  for  varying  heat  dissipation  characteristics  is  cur¬ 
rently  under  investigation. 

5.  Conclusion 

In  this  paper,  a  temperature-sensor-based  MPPT  technique  has 
been  proposed  for  TEG  energy  harvesting  systems.  The  presented 
system  operates  based  on  the  reference  voltage  for  MPP  generated 
by  the  feedforward  TEG  OCV  model.  The  proposed  MPPT  scheme 
can  track  the  MPP  of  any  TEG  system  without  current  measure¬ 
ment,  power  perturbation  or  source  disconnection,  which  enables 
optimal  energy  harvesting  at  any  temperature  or  load  condition 
with  a  much  simpler  and  much  less  expensive  circuitry.  The  MPPT 
scheme  can  be  implemented  with  a  minimal  number  of  discrete 
parts  and  requires  no  complex  hardware  or  software.  Moreover, 
it  is  inherently  stable  because  the  system  has  no  instantaneous 
power  measurement  feedback  loop  and  can  readily  be  used  in  con¬ 
junction  with  any  type  of  power  converter.  The  proposed  scheme 
has  been  validated  analytically  and  experimentally,  and  demon¬ 
strated  successful  performance. 
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